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ANALYTIC SOLUTIONS FOR THE BURGERS EQUATION 
WITH SOURCE TERMS 


G. I. MONTECINOS 


Abstract. Analytic solutions for Burgers equations with source terms, possi¬ 
bly stiff, represent an important element to assess numerical schemes. Here we 
present a procedure, based on the characteristic technique to obtain analytic 
solutions for these equations with smooth initial conditions. 


1. Introduction 

Exact solutions of hyperbolic balance laws are very useful to assess the perfor¬ 
mance of numerical schemes. Non-linearity of partial differential equation, as well 
as, stiffness of source terms are desirable features to be recovered by numerical 
schemes. The analytic solutions presented here contain these elements. 

Exact solutions for Burgers equations are generally obtained by separation of 
variables |T] , regularization techniques [3] , expansion methods , to mention but a 
few. Here, the strategy to solve these equations is based on the characteristic curve 
method, see HIS]. The solution is obtained in two steps. First, the solution of 
an Ordinary Differential Equation (ODE) called here, equivalent ODE, is obtained. 
This ODE is constructed by following the conventional characteristic curve method 
and contains the influence of the source term. Second, the solution of an ODE, 
called here characteristic ODE is obtained. This equation is defined in the x — t 
plane and the solution of the equivalent ODE is included. In this way the influence 
of the source term is present in the definition of the characteristic curve. 

To ensure the solvableness of the equivalent ODE, the existence of a primitive 
function for the reciprocal of the source term is required. Subsequently, if the 
solution of the equivalent ODE has a primitive function then the characteristic ODE 
is solvable. If these requirements are satisfied, analytic solutions can be obtained. 
However, a non-linear algebraic equation has to be solved in the general case, which 
require the initial condition to be a continuous function. 

This work is organized as follows. In section the procedure is presented. In 
sectionj^ the procedure is applied to the homogeneous Burgers equation. In section 
l^the solution for a linear source term is obtained. In sectionj^the analytic solution 
for a quadratic source term is obtained. In section the methodology is applied 
to solve Burgers’s equation with a more general non-linear source term. Finally, in 
section [3 the main results of this work are summarized. 


2. Analytic solutions for Burgers equations with a special class of 

SOURCE TERMS 


To start, let us consider a balance law in the form 


( 1 ) 


dtq{x,t) + doi 


\ _ 


= s{q[x,t)) 


q{x,Q) = ho{x) , 

where /io(x) is a continuous initial condition and s{q) is the source term which has 
to satisfy the properties of lemma 2.2 shown below. 
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( 2 ) 


df 


By following the characteristic curve method we define a curve x in the x — t 
plane, which satisfies 

‘'x(t) = q{x{t),t) , 
a;(0) = y , 

with y a constant value. This ODE will receive the name of characteristic ODE. 
On the other hand, we define 

(3) q{t) = q{x{t),t) ,h(0) := q{x{0),0) = ho{y) , 

where y is that given in With these definitions, 0 becomes an ODE given by 


(4) 


= s{q{t)) 


m = m. 

This ODE is called here, equivalent ODE. 

Now, let us see the following definition in order to call for ensuring the existence 
of solutions. 

Definition 2.1. The primitive of a function f{x) is a differentiable function F{x) 
such that 

(5) ^ ■ 

Next lemma sets the type of source terms which are considered in this work. 

Lemma 2.2. If s{q)~^ contains a primitive with respect to q, then the equivalent 
ODE is solvable. Additionally, there exists a function £{t,h(0)) such that 

§-^£{t,h{0)) = siq(t)) 

£i0,h{0)) = h{0) 

and q{t) = £{t, h{0)) . 

Proof. We integrate ([^ as follows 

fQit) ft 

(7) / s{q)~^dq= / dt , 

Jh(0) Jo 

as s{q)~^ has a primitive function, there exists G{q) such that 

(8) ^G(g) = s{q)-^ . 

Therefore, £{t, h(0)) is the solution to G(£) — G(h(0)) — t = 0 . So, if there exists 
the inverse function of G{q), which is denoted here by G~^{q), the function £ is 
explicitly given by 

(9) £(t,h(0)) = G-i(t + G(/i(0))) . 

In any case the exact solution is obtained as 


( 10 ) 


q{t) = £{t,h{0)) . 


□ 


Once the solution to the equivalent ODE is available and observing that h{0) := 
hQ{y), the characteristic ODE takes the form 

/l-|^ fMt) =£{t,ho{y)) , 

x{0)=y. 

The following lemma deals with the existence of solution for this ODE. 
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Lemma 2.3. If £(t, ho^y) has a primitive function IF(t,hQ(jj)), such that 
-^IF{t, holy)) = £{t,ho{y)) and IF{0, ho{y)) = 0 . Then, the characteristic ODE 
has the exact solution 

( 12 ) X = y + E(t,ho{y)) . 

Proof. Integrating the ODE ( [TT| ), we obtain 

(13) Jy dx = J^£{t,ho{y))dt. 

So, by using the properties of P{t, hQ{y)), the result holds. □ 


Remark 1. The value y, is a constant for the characteristic ODE {11). However, 
if values x and t are set in (12), there exist a constant y satisfying (12). Therefore 
we can identify a such constant by y = y{x, t) . 

Proposition 2.4. If s{q)~^ and £{t, ho{y)) have their respective primitive func¬ 
tions. The problem 0 has the exact solution 

(14) q{x,t) = £{t,ho{y)) , 
where y satisfies 

(15) X = yE{t,hQ{y)) , 

with IF(t,hQ{y)) the primitive of £(t,hQ{y)) with respect to t . 


Proof. The construction of this function is given by lemmas (2.21 and (2.3). Now, 
we are going to probe that q{x,t) solves ([^. 

By the chain rule 


d.q = • 


(16) 

dy^ dx 

Then 

(17) + + 

On the other hand, we have 

0 

1 = 


(18) 


— IIm. 4- JLt 
~ dt at 


dx 




Therefore, as is the primitive of £ and from (14l, we have 

-1 


(19) 




dt 


— 


‘ dx 


dho ^ 


i-i:F + q) = o. 


Finally, we note that 

(20) §^£it, ho{y)) = s{q{x, t)) 

and so, the result holds. 


□ 


Remark 2. Note that in the general case, (15) is an algebraic equation and the 
bisection method is used to solve them. The regularity requirements of the bisection 
method is that, functions have to be at least continuous ones, which is ensured by 
taking hQ{x) a continuous function. 

In the following sections we will present some applications of the strategy describe 
above. 
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3. Homogeneous Burgers equations 
In this section we apply the methodology seen in section to solve the problem 


( 21 ) 


dtq{x, t) + d^{ I = 0 


q{x,0) = ho{x) . 

Then, the equivalent ODE has the form 

= 0 , 


( 22 ) 


dt 

m = m ■ 


This ODE is solvable and the exact solution has the form q{t) = £(t, h{0)) = h(0) . 
Therefore, the characteristic ODE takes the form 


(23) 


f^x{t) = hoiy) , 
a;(0) = y , 


as £{t, h{0)) has the primitive function iF{t, h(0)) = th{0) , we obtain 


(24) 


x = y + tho{y) . 


However, ho{y) = ho(x(0)) = q(x(0),0) = q(x,t)- Therefore, (24) becomes 
(25) X = y + tq{x, t) . 


Hence the exact solution to (21), is computed as 
(26) q{x,t) = hQ{x-tq{x,t)) , 

which is the conventional solution of the homogeneous Burgers equation for general 
smooth initial conditions. 


Example 1. Let us consider the problem l[21\ l with Lq^x) = x . A simple inspection 
shows that, u{x,t) = is the exact solution. Following the present methodology 
we are going to obtain this exact solution. 

To start, let us note that equivalent ODE has the exact solution 

( 27 ) m = HQ) ■ 

Therefore, the exact solution to Ij21\ ) is given by 

(28) q{x,t) = ho{y) , 

with y satisfying the relationship obtained from the characteristic ODE, which by 
using the form of hg [y) has the form 

(29) x = y + thQ{y) = y{l + t) . 

So, by solving for y the sought solution is obtained 

(30) q{x, t) = ho{y) = y= . 

4. Burgers’s equation with a linear-source term 
Now, let us consider the PDE given by 

dtq{x, t) + da; = I3q{x, t) , 

q{x,0) = ho{x) . 

Then, the equivalent ODE has the form 


( 32 ) 


m = HQ) ■ 
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This ODE is solvable and the exact solution has the form q{t) = £{t,h{0)) = 
h{0)exp{/3t) . Therefore, the characteristic ODE takes the form 


(33) 


^^x(t) = ho{y)e^* 
x(0) = y , 


as £{t, h{0)) has the primitive function E{t, h{0)) = (^-g-^)h(O) , we obtain 

- 1. 


(34) 


x = y + 


-hoiy) ■ 


Therefore, the exact solution to lH) is given by 
(35) q{x,t) = ho{y)e^\ 


with y satisfying (34). 


Example 2. The function u{x,t) = is the exact solution of (31) with 

hQ{x) = , with (3 . 

To find this solution by using the present methodology we note that the equivalent 
ODE has the exact solution 


(36) 


q{t) = /i(0)e^‘ . 


On the other hand, the characteristic ODE due to the form of holy), has the exact 
solution 


(37) 


oht _ 1 


o/3t-l 


x = y + ho{y )—-— = y{l + ' 

Therefore, the exact solution to l[31\ l is given by 
(38) q{x,t) = ho{y)e^* , 


withy satisfying (37), which finally provides 
(39) 


q{x,t) = ho{y)e^* = 




- 2 + ■ 


Example 3. Let us consider the problem (31 ) in the interval [0,1] , with the initial 
condition ho^x) , given by 


(40) 


ho{x) = < 



2 

, X < 0.1 , 

2 + 

d 

1 o 

, 0.1 < T < 0.2 , 


3 

, 0.2 < T < 0.4 , 

3- 

(x—0.4) 

0.2 

, 0.4 < a; < 0.6 , 


2 

, X > 0.6 . 


In that case, the equation (34) is not explicitly solvable for y. In this work we use 
the bisection method to obtain y. 

shows initial condition (fO) and the solution of {31) at time tout = 0.2 

for (3 = —2 . 

5. Burgers’s equation with a quadratic source term 
Let us consider the partial differential equation 


dtq{x,f) + dj: 


qD,tf \ _ 


= I3qix,ty 


q{x,t) = ho{x) . 


(41) 
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— - q(x,0) 



1 ; 
0 


0.5 

X (space) 


1 




Figure 1. Burgers’s equation with linear source term. Initial con¬ 
dition (dash line) and the solution at time tout = 0.2 (full line) for 
/3 = -2 . 


The equivalent ODE has the form 

(42) 

which is solvable and the exact solution is 

(43) 


^ = m? 
m = h{o ), 


g(t) = £:(t,M0))= i_^j°^(0) . 


On the other hand, we note that the characteristic ODE 

f^x{t) = q{x,t) , 


(44) 


x(0) = y , 


for q{x,t) satisfying (431 has the solution 

ln{l - Ptho{y)) 


(45) 


x = y 


P 


Therefore, the solution of ( |41[ ) is given by 

ho{y) 


(46) 


qix,t) = 


1 - pthoiy) ’ 


with y satisfying (45). 


Example 4. Let us consider in the interval [0,1] , with initial condition 
(47) ho{x) = sin(27ra;) . 


With this choice, y is not explicitly obtained from (4-5). Therefore we use the 
bisection method. Eigure^ shows the initial condition (dash line ) and the solution 
at tout = 0.15 and for /3 = —2 . 
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Figure 2. Burgers’s equation with quadratic source term. Initial 
condition (dash line) and the exact solution at time tout = 0.15 
and for (3 = —2 . 


6. Burgers’s equation with a non-linear source term 
Let us consider the partial differential equation 

dtqix, t) + d, 

q{x,t) = ho{x) . 

So, the equivalent ODE 


(49) 

has the exact solution 


^<?(*) _ p3q{t) 
dt ^ 

m = h{o), 


(50) 


q{t) = £{t,h{0)) = - 


ln{< 


-ph(o) _ 


Pt) 


On the other hand, the characteristic ODE 

(51) - fS ’ ^ 

a:(0) = y . j 

with q{t) = q{x{t),t) satisfying ( [sl] ) has the exact solution 

(^e-fdhoiv) _ j3t)ln{e-^^^^y'i - pt) + pt + phoiy)e-l^^o^y'> 


(52) x = y + 


P^ 


Therefore, the solution to (48) is given by 

;„(e-/3'io(y) _ pt) 


(53) 


q{x,t) = -- 


P 


with y a solution of (52). 
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Figure 3. Bueger’s equation with exponential source term. Ini¬ 
tial condition (dash line ) and solution at tout = 0.26 (full line). 
Parameters /3 = — 1 and e = 10“^ . 


Example 5. Let us consider 


in the interval [0,1], with the initial condition 


(54) 

with w(x) = 


ho{x) = 2 
0.3 — X 


1 — w{x) 
2 

. Figure 


+ 


1 -|- w{x) 


(0.3 - xY +£ _ 

and the respective solution at time tout = 0.26 for (3 = —1. 


3 shows the initial condition for e = 10 


-4 


7. Conclusions 

In this work we have presented a procedure to obtain analytic solutions to the 
Burgers equation with a family of source terms aimed to provide a set of tests 
which is suitable for the assessment of numerical schemes. The family of source 
terms is that formed by functions containing primitive functions with respect to its 
arguments. By following the characteristic method, we have defined two ODE’s, 
which were called equivalent ODE and characteristic ODE. This family of source 
terms allows us the solvableness of these ODE’s and as a consequence the sought 
solution of the original PDE is found. In the solution of the characteristic ODE 
an algebraic equation, which depends on the initial condition, is obtained. The 
requirement of continuity for the initial condition allows us to find the root of this 
algebraic equation through the bisection method. 
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